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We present a new analytic approach for the study of late time evolution of linear test-fields, 
propagating on the exterior of black holes. This method provides a calculation scheme applicable 
to Kerr black holes (for which case no analytic calculation of the late time tails has been presented 
so far). In this paper we develop the new technique and apply it to the case of massless scalar 
waves evolving on the background geometry of a static spherically-symmetric thin shell with a 
Schwarzschild exterior. The late time behavior of the scalar field at null infinity is calculated, and is 
explicitly related to the form of (quite arbitrary) initial data. This reproduces the well-known late 
time power-law decaying tails. In an accompanying paper we apply our approach to the complete 
Schwarzschild black hole geometry, where we obtain the familiar inverse-power late time tails at null 
infinity, as well as at time-like infinity and along the event horizon. A detailed calculation of the 
late time power-law tails in the Kerr geometry, based on the same approach, will be presented in a 
forthcoming paper. 
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I. INTRODUCTION 

The gravitational field outside black holes created dur- 
ing a generic gravitational collapse relaxes to the sta- 
tionary Kerr-Newman field. This result, referred to as 
the "no-hair theorem" (see, for example, H) implies that 
outgoing radiation should carry away all the initial char- 
acteristics of the collapsing object, except for its mass, 
charge, and angular momentum. These three quantities, 
characterizing the Kerr-Newman field, are conserved by 
conservation laws; any other quantity will vanish by the 
time the gravitation field settles down on its stationary 
state. 

The underlying mechanism for this relaxation process 
was first demonstrated by Price || for the case of a nearly 
spherical collapse. Price analyzed the dynamics of mass- 
less integer-spin test fields, evolving on the background of 
a Schwarzschild black hole (SBH), and showed that when 
viewed from a fixed location outside the black hole, the 
waves die off at late time with an inverse power-law tail 
(in the Schwarzschild time t) , whose power index depends 
only on the multipole number I of the mode under consid- 
eration. Later it was demonstrated by analytic and 
numerical methods that the nearly spherical collapse ex- 
hibit late time decay tails also at future null infinity and 
along the event horizon. The formation of these tails 
was explained as due to back-scattering of the outgoing 
radiation off spacetime curvature at very large distances. 

The existence of late time tails was demonstrated for 
linear perturbations of both SBH and Reissner- 

Nordstrom |j) exteriors. Remarkably, numerical analysis 
of the fully non-linear dynamics of fields yields the same 
decay rates at late time, as for the minimally-coupled 
(linear) fields This, of course, encourages the ap- 



plication of the perturbative approach, even though the 
problem in discussion is non-linear in its nature. 

Ching et al. |J explored a wide class of asymptotically- 
flat spherically symmetric spacetimes, represented by ef- 
fective curvature potentials of the form (In r) 13 /r a (where 
j3 = 0, 1 and a > 2 are parameters). This class includes 
the SBH and Reissner-Nordstrom geometries. They 
found that, generically, the wave behavior at late time 
is not characterized by a strict power-law tail, but rather 
it has the form (In t)@x (inverse-power in t). In that re- 
spect, the Schwarzschild and Reissner-Nordstrom geome- 
tries represent a special subgroup of spacetimes. For the 
monopole moment of the scalar radiation, it was shown 
by Gomez et al. that the form of the tails (whether 
"logarithmic" or not) depends on whether or not the 
Newman-Penrose constant for the field vanishes. 

Recently [||, Brady et al. studied scalar waves 
dynamics in the non-asymptotically-flat exteriors of 
Schwarzschild-de Sitter and Reissner-Nordstrom-de Sit- 
ter black holes. Contrary to the asymptotically-flat ge- 
ometries, no power-law tails were detected in these cases. 
Instead, the waves were found to decay exponentially at 
late time. 

The aforementioned analysis by Ching et al. follows a 
technical scheme, first introduced by Leaver H, in which 
the linear waves are first Fourier-decomposed, then eval- 
uated in the complex frequency plane (see also 0). In 
this technique, the late-time tails are explained in terms 
of a branch cut in the Green's function in the frequency 
domain. The fact that the branch cut is due to the form 
of the potential at asymptotically large radius (see || 
for details) implies, again, that the tails originate from 
scattering off the curvature potential at large distances. 
This observation, in turn, suggests that the development 
of tails is independent of the existence of an event hori- 
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zon. Thus the tail phenomenon appears to be of a more 
universal nature: it may characterize the realistic stel- 
lar dynamics. This, indeed, was suggested by Gundlach 
et al. J3|,|| and Ching et al. §. In B, for example, the 
purely spherical collapse of a self-gravitating minimally 
coupled scalar field was studied. It was demonstrated 
numerically that in this case late time tails are formed 
even when the collapse fails to create a black hole. We 
shall further discuss this issue in the present paper. 

Historically, the study of wave dynamics outside black 
holes was motivated by the will to construct a detailed 
description of the relaxation process leading to the sta- 
tionary "no hair" state (For example, there was a need 
to verify that the event horizon is indeed stable under 
generic perturbations). Many recent studies, however, 
are oriented by the prospects of directly observing grav- 
itational radiation from astrophysical systems (e.g. with 
the LIGO observatory, see 111]]). Although the study of 
late time behavior of fields outside black holes has prob- 
ably no immediate observational implications (late time 
outflux from a realistic coalescing binary system, for ex- 
ample, would be orders of magnitudes weaker than the 
short pulse of radiation expected during the last few sec- 
onds of the merger) , it still addresses several important 
questions. For example, this study has crucial relevance 
to the exploration of the internal structure of black holes, 
as it provides the input for the internal wave evolution 
problem. Particularly, it has been argued Jljj that the 
form of the slowly decaying wave tail along the event 
horizon affects the strength of the mass-inflation singu- 
larity at the Cauchy horizon inside charged and rotating 
black holes. 

The study of waves evolving in curved geometries is 
also interesting on its own right, from a theoretical point 
of view: such waves, massless just as well as massive, 
do not propagate along light cones solely; rather they 
also spread inside them. Regardless of the presence or 
absence of an event horizon, it is this feature of the evo- 
lution which is responsible for the phenomenon of late 
time decay tails in curved spacetimes. 

In virtue of previous studies, we now have the follow- 
ing schematic picture regarding the dynamics of waves 
outside a nearly-spherical collapsing object: Consider a 
perturbation in the form of a compact pulse of radia- 
tion (gravitational or electromagnetic), somewhere out- 
side the collapsing object, emitted at some time during 
the collapse. This pulse may represent radiation emerg- 
ing from the surface of a collapsing object, as well as 
any other form of perturbation on the background geom- 
etry. We shall refer to it as the "initial pulse" . A static 
observer outside the black hole will then indicate three 
successive stages of the wave evolution. First, the exact 
shape of the waves front depends on the detailed form of 
the initial pulse. This stage (lasting a period of time com- 
parable to the duration of the initial pulse) is followed by 
a "quasinormal ringing" (QNR) relaxation stage, during 
which the waves undergo exponentially-decaying oscilla- 
tions with (complex) frequencies completely determined 



by the mass and electric charge of the central object p3| . 
Finally, as the QNR dies off exponentially in time, it 
leaves behind an inverse power-law decaying tail of radi- 
ation. During the last two stages, the details of the initial 
pulse affect the shape of the waves only through a global 
amplitude factor, hence the evolution during these stages 
is purely characteristic of the background geometry. 

Now, a realistic black hole formed by a generic grav- 
itational collapse is expected to spin, as do astrophys- 
ical stars (On the contrary, models of black holes with 
electric charge and/or cosmological constant are mainly 
hypothetical). Therefore it is natural to ask how the 
tail phenomenon is affected by the presence of angular 
momentum in the background geometry. One may sus- 
pect that the late time behavior of waves, evolving out- 
side a rotating Kerr black hole, would be qualitatively of 
the same nature as in the spherically-symmetric space- 
times. This is because the Kerr geometry asymptotically 
approaches that of a SBH at very large radius, and pre- 
sumably it is this far region of spacetime whose structure 
determines the form of the late time radiation. 

Yet, in spite of the obvious interest, no analytic scheme 
has been proposed thus far to describe wave dynamics 
outside rotating black holes. The basic obstacle is, of 
course, the fact that the axially-symmetric Kerr black 
hole possesses three non-trivial dimensions, instead of 
only two in the spherical cases. This makes both ana- 
lytic and numerical investigation significantly more com- 
plicated. Recently, there was an initial progress, with 
the introduction of a full (2 + 1 dimensions) numerical 
analysis by Krivan et al. firjjla], though much effort is 
still required is this direction. 

The new calculation approach 

We have developed an analytic calculation scheme en- 
abling one to analyze late time wave dynamics outside 
Kerr black holes. This analysis shall be presented in fjq] 
(see also (jl7|). The major goal of this paper, together 
with the accompanying paper (to be referred to as Pa- 
per II), is to demonstrate the applicability of our tech- 
nique in a simpler model, namely the evolution of scalar 
waves in the SBH exterior. This will serve several pur- 
poses: First, we shall be able to test our scheme against 
the firm results already obtained for this case by previ- 
ous studies. Secondly, many parts of the formalism to be 
developed shall be later directly employed in the analysis 
in Kerr. Finally, our analysis in Schwarzschild proves to 
be valuable on its own right, providing, in some respects, 
a more complete picture of the late-time wave behavior 
than already available. 

Basically, our analysis is composed of two major steps. 
In the first and more crucial one, a characteristic initial 
value evolution problem for the scalar waves is treated 
analytically, resulting in the construction of a late-time 
solution for the wave at null infinity. This calculation in- 
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volves the introduction of a special perturbative decom- 
position of the waves, followed by the application of the 
standard (time domain) Green's function technique. In 
the second step we then use a simple late time expansion 
of the wave near time-like infinity, in order to obtain the 
late time behavior of the wave at any constant radius, 
including along the event horizon. This second step is 
made possible only after the waves' form at null-infinity 
is derived by the first step of the analysis. 

To introduce the analysis at null infinity (i. e. the first 
step mentioned above) in a clear and more instructive 
way, we incorporate in this paper a simple toy-model, in 
which the scalar waves are taken to evolve in the grav- 
itational field induced by a spherically-symmetric thin 
shell of matter. This configuration possesses the same 
spacetime structure at large distances (outside the shell) 
as in a complete SBH manifold of the same mass, while 
its small-r structure is much simpler. That would re- 
duce the amount of technical details to deal with when 
developing our calculation scheme, and may also enable 
one to push the analytic calculation to a further extent, 
while leaving the essential features of the analysis un- 
affected. Indeed, we later show in paper II, consider- 
ing the complete SBH manifold, that the asymptotic late 
time behavior calculated at null-infinity is unaffected by 
the structure of spacetime at small radii, thus it is essen- 
tially the same in the complete SBH manifold as it is in 
the shell toy-model. 



Arrangement of this paper 

In this paper we study the thin-shell toy model. To 
that end we introduce our calculation method, which we 
call "the iterative scheme" , and apply it to calculate the 
behavior of the wave at null-infinity in this model. 

The paper is arranged as follows: In sec. || we for- 
mulate the problem of wave evolution in the shell model 
as a characteristic initial-value problem. The iterative 
scheme is p resen ted in sec. III. It is applied in sections 
[i"v| through VIII , revealing the power-law pattern of the 
late time wave decay at null infinity. The amplitude of 
the wave at late time is explicitly calculated, expressed 
in terms of a (quite arbitrary) initial data function. 



II. INITIAL- VALUE FORMULATION OF THE 
WAVE EVOLUTION PROBLEM 

A. The shell model 



ds 2 = - f{r)di 2 + f^i^dr 2 + r 2 {d9 2 + sin 2 9d<p 2 ), (1) 

where i, r, 9 and tp are the standard Schwarzschild co- 
ordinates and /(r) = (1 — 2M/r). The region inside the 
shell, r < R, is flat: 

ds 2 = -dT 2 + dr 2 + r 2 (d0 2 + sin 2 9dtp 2 ). (2) 

Expressed in these coordinates, the geometry suffers a 
discontinuity at r = R (the metric functions g rr and gu 
jump through this surface. Also the time coordinates do 
not agree on both sides of the shell). 

To allow a continuous representation of spacetime ge- 
ometry, we define the coordinates 



f(R)- 1/2 r 



,r<R 



and 



R + Rf(R)-^ + 2 M\n(j r 3Ml) ^>R 



f{R)- 1 ' 2 T ,r<R 



t 



,r>R. 



(3) 



(4) 



In terms of the coordinates i*, r + , 9 and (f> the geome- 
try is described as a single continuous manifold. (These 
coordinates form an example of a 'natural' set of coordi- 
nates through the 'surface layer' r — R, as discussed by 
Israel in The jump in the derivatives of the metric 
functions through the shell is then related to the surface 
energy tensor of this layer.) Note that outside the shell, 
the radial coordinate r* is the usual Schwarzschild's 'tor- 
toise' coordinate, satisfying dr^/dr = / -1 (r). 



B. The initial- value problem 

We consider the evolution of initial data, represent- 
ing a generic pulse of massless scalar radiation, on the 
fixed background of the shell described above. The scalar 
field is assumed to satisfy the (minimally-coupled) Klein- 
Gordon equation^] 



0, 



(5) 



where $ represents the scalar wave. The structure of 
spacetime affects the evolution of the scalar wave through 
the covariant derivatives, denoted in Eq. (j^) by semi- 
colons. 

Decomposing the field in spherical harmonics, 



We consider a static spherically-symmetric thin shell 
of matter, of a mass M and some radius R > 2M. The 
parameter R should be regarded as being of order ~ 2M 
(say, R — 3Af). For this configuration, the exterior vac- 
uum region r > R is a part of the Schwarzschild space- 
time, described by the line clement 



*We comment, however, that other wave equations are also 
possible. For example, to assure conformal invariance, one 
should rather use the equation □ $ + |^ ( J > = 0, where R 
is the Ricci scalar. This equation reduces to Eq. (||) when 
considering waves in vacuum. 
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*(<,r,0,¥>)=x; E ^(*.»-)^m(^v), 



(6) 



and substituting in the wave equation (|^), we obtain an 
independent equation for each of the components <f> (t, r). 
(We use the superscript I to denote the multipole number 
of the mode under consideration.) 

A convenient form for the wave equation may 
be obtained in terms of a new wave function 

(t, r) = r(/) 1 (t, r) . To that end we introduce the double- 
null (Eddington-Finkelstein) coordinates v = t + and 
u = t — r* . The equation governing the evolution of the 
l-mode of the scalar wave then reads 



Kv + V l (r)* l = 0, 



in which the function 

Hi- 



V l (r) 



IK 2 



+ 



2M 



,r <R 
,r> R 



(7) 



(8) 



serves as an effective potential for the scalar field. The 
form of the potential function completely determines 
the effect of spacetime curvature on the wave evolution 
[terms proportional to M], as well as the centrifugal ef- 
fect on the non-spherical (I > 0) modes of the wave [terms 
proportional to 1(1 + 1)]. 

To set-up the initial value problem for the scalar ra- 
diation, one should supply initial data for each of its 
modes. Since in spherical symmetry each mode evolves 
separately, it would be sufficient to analyze the evolution 
of a single /-mode of the wave from a given initial data. 
We shall specify this initial data on two characteristic 
(null) surfaces, as illustrated in figure [T|J Specifically, 
we shall consider initial data in the form of some com- 
pact outgoing pulse, specified on the ingoing null surface 
v = 0|] 



0) = 



= 
T l (u) 



(9) 



where T (u) is some function of a compact support be- 
tween retarded times u = u$ and u = u\ (for brevity 
we henceforth usually suppress the /-dependence of the 
functions under consideration). We shall assume that 
| till > namely that the support of the initial func- 
tion T(u) is completely outside the shell. 

The choice of compact ('localized') initial data proves 
to be convenient for the purpose of calculation, ft also 



T For a characteristic initial value problem to be well posed, 
it suffices to specify only the value of the field on the initial 
characteristic surfaces (and not its derivatives). 

*The choice v — does not limit the generality of the ini- 
tial setup, due to the time translation invariance of the back- 
ground geometry. 



becomes useful when we later try to characterize the way 
in which the late time behavior of the wave depends on 
the location of the initial pulse. For that reason, it is 
instructive to consider the case |iio| 3> u\ — uq, for which 
Mo becomes a single parameter describing the location of 
the initial pulse. 




FIG. 1. The set-up of initial data. Shown is the confor- 
mal diagram representing spacetime geometry in the shell model. 
The spacetime region outside the shell (r > R) is a part of the 
Schwarzschild manifold, while the interior of the shell (r < R) is 
flat, u = «o an d v = are two initial characteristic (null) surfaces. 
An initially compact outgoing pulse is specified on v = 0. 

To complete the specification of the initial value prob- 
lem, a boundary condition should be set at the origin of 
coordinates. From a physical point of view, a wave start- 
ing out of regular initial data and traveling outside the 
shell, should not develop any singular behavior. Also, 
the field should remain regular (in fact, continuous) at 
the shell itself, since the jump in the potential function 
through this surface is finite in magnitude . Then, no 
singularity is expected to occur during the free evolution 
of the scalar field in the Minkowski interior of the shell. 
Enforcing regularity on cj> then automatically implies that 
the revised wave function ^ = r<f> should vanish at the 
origin r* = r = Q. To rule out unphysical divergent so- 
lutions, we hence impose the supplementary boundary 
condition 



*(n = 0) = 0. 



(10) 



In appendix |A| we show that the evolution equation (Q), 
together with the initial conditions (|^) and the bound- 
ary condition (pTj|), establish a well defined characteristic 
initial value problem for the scalar field, with a unique 
solution anywhere inside the light cone of the initial data. 



III. THE ITERATIVE EXPANSION 

Our goal is to calculate the late time behavior of \& 
at null-infinity (namely, for v — * oo, u 3> M). In what 
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follows we introduce an analytic scheme to construct this 
solution. This technique is based on the iterative expan- 
sion, to be defined now. 
First, define 



IV. DERIVATION OF * 



o 



Vo(r) 



1(1 + 1) 
4 r 2 



and 



SV(r) = V{r)-V {r), 



(11) 



(12) 



with the potential V(r) given by Eq. (||). 

Then, consider a decomposition of the wave function, 

* = f + *i + *2 + ---, (13) 
such that the components jy obey the recursion formula 



N 



,N = 

-(SV^n-! , A > 



(14) 



(15) 



and satisfy the initial conditions 

q N (u = u ) = (VA > 0) 

T , / ^ f r(u) ,n = o 

and the boundary conditions 

*w(n=0) = (VA>0). (16) 

Formally summing Eqs. (14jl^Jl^) over N, we recover 
Eqs. fl^,P,|lO|) for the complete wave "P. This suggests 
that if the sum ([b|) converges, it should yield the correct 
function ^. 

Equations (|ljjl^,|l^) constitute an infinite hierarchy 
of initial value problems for the functions ^n- Each of 
the evolution equations consists of a simple homogeneous 
part, which is the free evolution equation in Minkowski's 
spacetime. In addition, each of the functions ^>n (ex- 
cluding ^o) has a source term proportional to the pre- 
vious function in the series. Therefore, in principle, if 
the solution for \to is found, and the appropriate (time 
domain) Green's function is constructed, then we should 
be able to solve for the functions one by one, using 
the standard Green's function method. 

Analytic arguments, as well as numerical simulations 
(both to be presented later) strongly suggest that at null- 
infinity the series ([l3|) does converge to the correct func- 
tion ^. Moreover, in the case \uq\ ^ M, which will con- 
cern us here, we find (both numerically and analytically) 
that the late-time behavior of \1/ at null-infinity is well ap- 
proximated by (the corrections coming from A > 2 
are smaller by a factor proportional to M/uq). Thus, 
in our scheme, understanding the behavior alone will 
suffice to determine the essential features of the late-time 
dynamics at null-infinity. 

In the sequel we derive an exact analytic expression for 
\l/o and for the Green's function , and use these results to 
calculate the late time behavior of at null-infinity. We 
then discuss the contributions coming from higher orders 
of the iterative expansion. 



By definition, admits the homogeneous wave equa- 
tion 



(17) 



where Vq is the purely centrifugal potential defined in 
Eq. (|TT|) . This equation is supplemented by the initial 
conditions 



, U = Uq 

T(u) ,v = 



(18) 



and by the boundary condition specified in Eq. ([16|). (re- 
call that r is a function of compact support, representing 
the initial pulse of scalar radiation.) 

In fact, is nothing but the solution of the analogous 
wave evolution problem in Minkowski spacetime||. To see 
that, notice that inside the shell Eq. ( |l7j ) is invariant un- 
der the transformation from the , t* coordinates to the 
usual flat-space coordinates r and t. Outside the shell, 
however, only the functional form of ^o(u, v) is the same 
as in Minkowski, while the dependence of the character- 
istic coordinates u, v on the flat space coordinates r, t 
differs. 

In appendix ^ we show that a solution to Eq. ([l7]) , 
subject to both the initial conditions ( |l8| ) and the bound- 
ary condition at r = 0, is unique. This solution reads 



i 

toS<«<0)^4l r 



9o ( u ) 



- u 



® (u > 0) = 0, 
in which the coefficients A l n are given by 



(2l-n)\ 
n\(l — n)\ 



(19a) 



(19b) 



(20) 



The function go(u) (with its parenthetical supcrindiccs 
indicating the number of times this function is differenti- 
ated) is, by Eq. (|l8|), the solution to the inhomogeneous 
ordinary equation of order /, 



71=0 



.4 



1 au ^ ' = r(u) 



(21) 



subject to the initial conditions g^\u Q ) = for all < 
n < I — 1. (Due to the compactness of F(it), we then 
automatically have also g^(uo) = 0. Thus vanishes 
along the ray u — uq, as necessary.) 



* In this respect, the iterative decomposition is actually "an 
expansion of Schwarzschild about Minkowski" . 
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The solution for go is given by 
1 f u / u \ '+i 



go(u < 0) 



(/-!)! 



(ni \ b-\-l 
— ) (u - u'y^Tfu^du' 
u' J 



(22) 



for I > 1, and simply go(u) — T(u) for / = OQ. It is easy 
to confirm this result by a direct substitution, noticing 
that a solution to the homogeneous equation correspond- 
ing to Eq. ( f2l|) has the general form 



g (ui <u<0) = 



,1 = 



J2 l k=iiku l+k ,i>o ' 



(23) 



where are constant coefficients. Note that this is also 
the form of go at u\ < u < (for which retarded times 
Eq. (^lj) becomes homogeneous) , with the coefficients 
being certain functionals of the initial data function T, 
that can be easily constructed by comparing Eqs. ( p^ ) 
and@. 

We find that during retarded times uq < u < tti, at 
which the initial pulse is "turn on" on the initial ingo- 
ing null surface, the behavior of ^o(u) depends on the 
detailed structure of the pulse. At later (yet negative) 
retarded times, u\ < u < 0, the function go(u) takes the 
s imp le form u l+1 x polynomial of order I — 1 in u. Eq. 
(19a) then implies that the wave itself dies off (as 
~ u, generically) towards retarded time u = 0. Excep- 
tional is the I = mode of the wave, which vanishes right 
after the initial pulse 'ceases' on the initial ray (that is, 
it vanishes identically at u > u\). 

In the later zone of spacetime, u > 0, the only solution 
for satisfying both ^q(u = 0) = Ofj and ^a(r* 
i) : = 0, is the trivial null solution ( |l9tf ) [see figure |). 
The uniqueness of this solution is guaranteed in virtue of 
the discussion in Appendix [X| 

We can now unify both Eqs. ( |l9| ) by defining 



g (u > 0) = 0. 



(24) 



It is also convenient to define go(u < 0) = 0, which makes 
go{u) a function of com pact support between u = uq and 
u = 0. Then, Eq. (19a) suffices to describes the wave 
anywhere. 



**In the scalar model, the monopole mode is radiative as 
any other mode. 

t ' Here we assert that the wave is continuous through the 
ray u = 0. For, as can be easily verified, a discontinuity along 
this ray must result is a violation of either the field equation 
or the boundary condition. The continuity of the wave is also 
reasonable from the physical point of view. 




FIG. 2. The evolution of \l/o (representing the scalar wave in 
Minkowski spacetime). For a compact initial pulse 'turned on' 
during retarded times uo < u < u\ < 0, the support of \PoM is 
confined to uq < u < (or uq < u < u\ for 1 = 0). The region in 
which has non-zero amplitude (for I > 1) is represented by the 
shadowed area in this conformal diagram. 

We conclude that for any mode of the scalar radiation, 
the evolution of is "cut off" not later then at u = 0. 
This somewhat surprising feature of the scattering off 
the purely centrifugal potential originates, in the 1+1 
representation, from a destructive interference between 
ingoing and outgoing wave fronts at the origin = 
(see appendix [Aj for details). This may be more eas- 
ily understood in the corresponding 3+1 picture, where 
compactness of the initial pulse directly leads to com- 
pactness of the wave fronts. The compact (in terms of 
retarded time) region in which survives is represented 
by the dark-colored area in the diagram of figure ^| 

Fi nally , we emphasize one further feature of ^o- By 
Eq. (19a), we have at null-infinity 



^(u) = g ( l \u) 



(25) 



in which 



(u) stands for ^!q(v — > oo). Now, the fact 
that go(u) is compact with respect to retarded time, 
means that any of its derivatives is compact as well. We 
thus find that the I integrals (*8°) (_1) , , 
carried out over all values of u at null infinity, all vanish: 



W) ( " n) = gt n) 



(l<n<0- 



(26) 



This feature shall appear to have an important impact on 
the form of the inverse power-law behavior of the scalar 
wave at late time. 

Obviously, since vanishes identically at u > 0, it 
does not contribute to the late time radiation. Rather, 
it serves as a source to higher terms in our iterative ex- 
pansion, as we show below. 
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V. CONSTRUCTION OF THE GREEN'S 
FUNCTION 

In order to calculate the next terms in the iterative 
expansion (that is to obtain solutions to the hierarchy of 
inhomogeneous equations (|l4|) for N > 1), we shall use 
the standard Green's function approach. To that end, we 
first need to obtain the Green's function corresponding 
to the operator d v d u + Vq. This is the purpose of this 
section. 

The (retarded) Green's function G(u,v;u' ,v') shall be 
defined as a solution of the equation 



G tUV + V^G = 5{u - u')6(v - v'), 



(27) 



subject to the causality condition G = outside the fu- 
ture light cone of the delta source point at («',«'). In 
the sequel it will become evident that this condition, to- 
gether with an appropriate boundary condition, specify 
a unique solution for G. 

To construct this solution we consider separately the 
two distinct regions of spacetime indicated in figure 
which arc defined with respect to given a source point 
(representing a source 2-sphere in 3+1 dimensions) at a 
certain location (u',v'). These are (I) the region inside 
the future light cone of (vf, v') not causally influenced 
by the origin (that is v > v' , vf < u < v'); and (II) 
the region inside the future light cone of (u',v') which is 
causally influenced by the origin (that is u > v'). 

First consider the Green's function in region I. We im- 
pose causality by writing 



(28) 



G = G9{u-u')6{v-v'), 



where G is a solution to the homogeneous equation (|2_ 
at u > vf and v > v', and 9 is the standard step function 



By analogy to Eq. (19a) we then have the solution 



G( U ,V-u',v') = Y. A n 



9 G \ u ) 
(v — u) l ~ n ' 



n=0 



(29) 



in which the coefficients A l n are those given by Eq. (20), 
and where the function ga{u;u' ,v') is yet to be deter- 
mined. (Alternatively, we could have expressed this so- 
lution in terms of a function ha(—v) instead of g(u). Of 
course, the unique solution for G, to be derived below, 
would have been the same). 

Inserting Eq. (pj|) into Eq. (|27j), we find that we must 
have 



G(u = u') = G(v = v') = 1. 



(30) 



With Eq. ( p9|), the equality G(v = v') = 1 takes the 
form of Eq. ( pi|) , in which we first make the replacements 
.9o ~ > <?Gj u ~± i u ~~ v ') and T(u) — > 1. Consequently, we 
find by Eq. ( ^2|) that (in region I) 

9G (u;u',v') = [ U ( 4^7 («" u") l ^du". 



(/-I)! 



(31) 



To calculate the integral we make use of the formula 



(x - XxY 



[x - x 2 ) n2 



dx = 



(na-i-2)! (x - zi)" 1 " 



(m - j)!(na - 1)! {x - na)" 2 ^'" 1 



(32) 



which is valid for any two natural numbers n\ and n 2 
satisfying n 2 > n\ + 2. (we shall need this integral sev- 
eral more times in the sequel). With some additional 
algebraic manipulations we can then finally obtain 



g G (u;u',v') 



(v' — u)(u — u') 



(v' - u>) 



(33) 



It is easy to verify that with this result, we have also 
G(u — u') — 1, as necessary. 



r=0 





FIG. 3. The Green's function for a given source 2-sphere at 
(u' ,v'). G differs from zero only in region I (the shadowed area) . 
It vanishes in region II. 

Next, consider region II, where each of the functions 
^jv is subject to the boundary condition ^n^* = 0) = 0. 
We impose this condition by requiring that the Green's 
function connecting any source point to the origin should 
vanish: 



G(r* = 0;u',v') = 0. 



(34) 



Note that this boundary condition fails to be valid at the 
single point u = v = v', where its value, dictated by the 
evolution in region I, is found by Eqs. ( |29| ) and (31) to 
be unity. Nevertheless, the boundary condition for the 
various functions "J at is still satisfied at all values of u, 
as will become clear later. 

In appendix [a| we show that, with the above bound- 
ary condition imposed, the only solution to Eq. (j27j) in 
region II is the null solution, namely G(u > v') = 0. We 
conclude that the Green's function is given by 



G(u, v; u , v ) 



£4 



9g ( 



u: u , v 



x 9{v 



2 (v-u) l - n 
v')9{u - u')9{v' 



u), (35) 
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where the function go is the one specified in Eq. d33|).P| 
In figure ^ we have indicated the region of spacetime 
to which the Green's function G(u,v; ,u',v') carries the 
influence of a given source point located at (u',v'). In 
practice, however, we shall be interested in calculating 
the functions ^n(u, v) (for N > 1) at a certain evaluation 
point (u,v). This will involve integration over all point 
sources (u', v'). The region in which sources influence the 
behavior of the wave at a given evaluation point (it, v) is 
indicated in figure ||. 



i+ 




FIG. 4. The Green's function for a given evaluation point at 
(u,v). The wave ^>jv at that point is influenced only by ^jv-i 
sources inside the shadowed area. 

VI. CALCULATION OF *i AT NULL INFINITY 

We have seen that the first element of our iterative 
series, namely *S>o, do not contribute to the overall late 
time radiation. Rather, it serves as a source to higher 
terms in the series. The first contribution to the late time 
radiation comes from ^i. In what follows we shall give 
a detailed calculation of Specifically, we will obtain 
an analytic expression for $1 at null- infinity, evaluated 
to the leading order in M/u. The result will be highly 
significant, since this term will appear to be the dominant 
constituent of the overall wave ^ at null infinity. We 
shall discuss this issue later, when we analyze the higher 
(N > 2) terms of the iterative expansion. 

For we have, by definition, 

*i,uv+V V 1 = -(5V)V , (36) 



** Note that the Green's function suffers a discontinuity 
along the outgoing ray u = v . This discontinuity, which is of 
order unity, originates at the delta source point (u',v') and 
travels along the ingoing ray v — v'. When this discontinuity 
encounters the origin, it is reflected back along u = v . Some 
subtleties related with the appearance of this discontinuity 
are discussed is appendix 



with the initial conditions ^>\(u = uq) = ^i(v = 0) = 
0, and the boundary condition ^(r* = 0) = 0. The 
functions Vq(t) and 8V{r) were given in Eqs. ( |Ti"| ) and 
(|l2l), respectively. 

The solution can be formally written as 

*i(tt, v) = - jdu'j dv' G(u, v; u', v')SV{u', v')V (u', v') 

(37) 

[with ^>o(v' < 0) = understood], where G is the Green's 
function in Minkowski's spacetime, given in Eq. (|35|). 
This form manifestly admits the above initial conditions. 
Also, the boundary condition is clearly satisfied, as nec- 
essary. Notice that although Eq. ([34]) fails to hold at 
u — v = v' , still the wave ^1 vanishes at that point [be- 
cause the v' integration in Eq. ( |37| ) has no support in this 
occasion.] 

The region of integration in Eq. (^) , for a given eval- 
uation point (it, v), is represented by the rectangle shown 
in figure || Sources outside this rectangle do not influ- 
ence the behavior of ^1 at (u, v). In addition, we have 
found previously that the support of ^ is confined to 
retarded times uo < u < 0. Hence, effective sources to 
^1 arc located only at the region indicated in figure || 
by the intersection of the two shadowed areas. We thus 
observe that when evaluated at late retarded time u, ^1 
is influenced only by sources at large radii. As u — > 00, 
it is only the asymptotically- far region of spacetime that 
affects the behavior of ^>\. 




FIG. 5. Region where ^0 sources affect at null infinity (in- 
dicated by the intersection of the two dark colored areas). If \ti is 
evaluated at late retarded time (u 3> M) then effective sources are 
confined to large distances (r' 2> M). 

We shall now evaluate Eq. (|37]), in order to yield an 
analytic expression for at null infinity, that is in the 
limit v — > 00. In this limit, the Green's function reduces 
to the form 

gl 

G(u,v -» oo;u,v) = -=-j[gG(u;u ,v )]. (38) 
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with only the n = I term surviving in Eq. ( p5[ ) , and where 
g G is given by Eq. @. 

We then notice that the /-derivative of go with respect 
to u can be taken out of the double integral in Eq. (37). 
The resulting surface terms all vanish, due to the factors 
(u — u') 1 and (v' — u) appearing in the numerator of go 
[see Eq. (Q)]. Therefore, in analogy to Eq. (25), we can 
now have 



(39) 



where we define 



It pi) 

gx(u,v) = - I du'l dv' g G (u;u' ,v')5V(u' ,v')$ (u' ,v'), 

J UQ J U 

(40) 

and where and g%° stands for the value of these func- 
tion at null infinity -> oo). It would hence be sufficient 
to calculate g^° in order to immediately obtain ^J . 

To allow explicit integration in Eq. (|40f ) , we must first 
express the function SV(r) in terms of the null coordi- 
nates. This cannot be done explicitly outside the shell, 
since the function r(r*) is implicit at r > R. However, 
in terms of the M/r* expansion we can write (Recalling 
that the parameter R is of order M) 



5V(r* > r R ) 



Mn2f» 
IP 







M 2 (lnf 



z ^2 



(41) 



where f * = r*/(R— 2M) , a is some constant (depending 
on R and M), and 



b = 8Ml(l + l). 



(42) 



Throughout the rest of this paper, the symbol "~" , when 
appearing over a character, shall always indicate the ratio 
of that quantity to R — 2M. 

Note that inside the shell we have SV = by defi- 
nition, and so the region r* < r*(R) = Rf(R)" 1 / 2 is 
automatically excluded from the domain of integration 
in Eq. @. 

The various terms in the M/r* expansion of SV con- 
tribute to "Pi in an additive way [via Eq. J40|)]. These 
contributions may be, in principle, calculated for each of 
the terms in separate. However, the calculation to follow 
suggests that to the leading order in M/u, it is only the 
leading order of SV in M/r* which contributes to at 
null infinity at late time . We shall therefore focus now 
on the contribution coming from this leading order. 

With the explicit form of g G and ^o, an d taking SV to 
the leading order in M/r* , Eq. (f40|) reads at null infinity 



■WA l n rd U 'g<?\u')(dr 

L n=0 Ju ° Ju 



,{v' -u) l {u-v!) 



(v' - u') 



x [a + bln(v' -u')]6[ri - n(R)}, 



(43) 



where the fact that SV = inside the shell accounts for 
the step function appearing in the integrand. 

To proceed, we discuss separately the cases I = and 
I > 1. First, consider the case I > 1. For these modes 
we show that the 'r~ 3 potential' source (i.e. the one pro- 
portional to a) contributes nothing to g^° at u > 0. An 
analytic expression for the contribution coming from the 
'logarithmic potential' source (the one proportional to b) 
shall then be derived. 

To show that contribution from scattering off the r~ 3 
potential vanishes (in the case I > 1), consider the term 
in d43| ) proportional to a. Integrating over w', with the 
help of Eq. (|32"|), this term yields 



rdu'(u-uT- 2 9 { o n) (u'), 

n=0 Ju « 



(44) 



with a vanishing contribution from the upper boundary 
of the ^'-integration, and where only the j = I term has 
survived in Eq. (p2l). The coefficients B l n are given by 



D 



l-n + 1 



n\(2l-n+l)(2l-n + 2)' 



(45) 



We now observe that in Eq. (^) all terms correspond- 
ing to 2 < n < I vanish independently for u > 0, since 
go(u) is compact. This can be verified by integrating each 
of the terms in the sum by parts n — 2 successive times 
with respect to u' , and then using Eq. ( |26| ) to find that all 
resulting surface terms vanish. Moreover, one also finds 
that the two remaining terms of Eq. (44), corresponding 
to 11 = and n = 1, add to zero. This is easily shown by 
noticing that B l Q = B[ for any I > 1. 

We conclude that all terms proportional to a in Eq. 
( p3| ) vanish (for u > 0), and thus that scattering off the 
r * potential does not affect the late time behavior of the 
I > 1 modes at null-infinity. This seemingly odd feature 
becomes clear by the following simple argument. Con- 
sider the solution \1/q to the wave equation ( |l7| ) , in which 

we take the potential to be Vq — j^^i- In addition, 
we take to satisfy the same initial conditions as ^0, 
and the boundary condition 5'5( r * = e) = 0. Clearly, we 
have ^(e — + 0) = ^o- Differentiating the wave equation 
for with respect to e and taking the limit e — >• 0, we 
obtain 



where "fg = lim e ^( 



V % 



9*o 



1(1 + 1) 
2rl 



(46) 



Comparison of Eq. ( f!6[ ) with 
the wave equation (|36| ) for ^1, shows that should be 
proportional to ^1 calculated with SV ~ r^ 3 . (Both 
functions, and ^>' , are subject to the same initial 
conditions: on the initial surfaces we have, by definition, 
^1 = and also 'I'q = 0, since the initial conditions 
chosen for ^ are e-independent.) Now, it is not difficult 
to verify that ^g vanishes identically at u > (as does 







\&o)0 Therefore the contribution to \&i due to scattering 
off the r~ 3 potential must vanish identically at u > 0, as 
demonstrated above for by explicit calculation. 

Next, we have to consider contributions to com- 
ing from scattering off the 'logarithmic' potential (i.e. 
contributions to the integral in Eq. (^) due to terms 
proportional to b) . The details of this calculation are left 
to appendix |^, where we obtain the simple result 

9?{u)=-2M [°toM d ^ (47) 

The case I — : The monopole case is especially 
simple to handle, as in this case no logarithmic terms are 
involved in the calculation (the coefficient b in Eq. ( ffl| ) 
vanishes), and the Green's function is simply unity. For 
I = we have no summation over n in Eq. ( |43| ) , and have 
a = 4M. This equation then reads 

rf» = -AM £ du> £ dv'^^ (for I = 0), 

(48) 

which directly leads to Eq. (|47|) . Therefore Eq. @) is 
correct for any of the modes. (Note the interesting result 
that g^° depends on I only through the explicit form of 

go)- 

Thus far we were approximating the potential SV by 
its large form, indicated in Eq. (pd]). In principle, to 
obtain the correct expression for gf, the complete po- 
tential SV should be considered. Each of the terms in 
the M/r* expansion of the potential SV contributes ad- 
ditively as a source to gi(u). These terms are all pro- 
portional to r~ k2 (lnf*) fcl , where k\ and k 2 are natural 
numbers, satisfying k 2 > 3 and < k\ < k% — 2. Above 
we have given a complete analysis concerning the leading 
contribution, namely the one corresponding to k 2 = 3. 
We can show, following a completely analogous analytic 
treatment, that the contribution from the above general 
term of the SV expansion (with k 2 > 3) is dominated at 
late retarded time by 

gT(u;k 1: k 2 )^M k ^ 2 f g^u' ^y^'X du'. (49) 

(Note that the k 2 — 3 case is special, in that in this 
case there is no logarithmic dependence involved in the 
resulting form of g^° ■) We thus conclude that omitting 
the contributions from terms other than the leading term 
in SV, has no effect on g%°(u) to the leading order in M/u 
and in uq/u. 



^ A straightforward calculation shows that *I/o dies off at 
u > 2e as ~ exp[— it/(2e)]. This implies the vanishing of 
d^o/de at u > in the limit e -> 0. 



Hence, in virtue of Eq. fl47| ) we obtain, to the leading 
order in M/u and in Uq/u, 

g^(u^>u ) = -2MI Q u- 2 , (50) 

in which we have defined 

/ g (u')du'= / g (u')du' = I . (51) 

Juq J— OO 

(The first equality is valid in view of the compactness of 

go(u).) 

By Eq. ( p9|) we can now finally obtain (to the leading 
order in M/u and in uq/u) 

^{u > u ) = 2(-l) l+1 (l + 1)\MI ir 1 - 2 . (52) 

Note that initial data information is manifested at null 
infinity at late time only through the single integral Iq. 
Also, note that in our approximation there is no reference 
whatsoever in the form of to the radius R of the shell. 

Both the power index of the radiation tail and the am- 
plitude coefficient deduced above, are the same as those 
obtained by Gundlach, Price and Pullin (GPP)(Eqs. (11) 
and (18) in ||) for the late time behavior of the 'com- 
plete' scalar wave at null-infinity, , using a different 
approach. Yet, our result regards only one element of 
the 'complete' wave, namely Vt^ . To reveal the overall 
behavior of the wave we must analyze each of the terms 
in the iterative scheme, then sum up their contributions. 
Thus, in order for our result to coincide with that of 
GPP, '5^° and the higher terms in our expansion 'must' 
be negligible at null-infinity. In the following chapter we 
argue that this, however, is not the case in general. We 
will show that when considering generic initial date, one 
finds additional unnegligible contributions to at late 
time. These contributions, coming from ^3, ■ ■ ■, shall 
have the form M 2 Uq 1 u~ l ~ 2 , M 3 Uq 2 u~ l ~ 2 , etc. In this 
respect, the result of GPP is not strictly correct: the 
amplitude calculated in H should fail to represent the 
amplitude of ^>°° for a generic initial pulse. However, as 
we indicate later, the result of GPP, as well as the form 
of calculated above, approximates the 'real' behavior 
of the scalar wave in the case where the initial compact 
pulse is confined to large radii, away from the highly 
curved region of spacetime. Mathematically speaking, it 
will be argued that the equality ^1 = ^ holds at null- 
infinity to only the leading order in M/uq. 

It is interesting to mention that other previous at- 
tempts to calculate the late-time behavior of scalar waves 
in Schwarzschild |^,|6|,p|jic|l , all yielded an expression pro- 
portional to the mass M, omitting possible contributions 
from higher powers of M. To our best knowledge, no 
analytic method has been proposed to enable calculation 
of the late time radiation resulting from initial date lo- 
cated at small radii, where curvature is large. In this 
respect, our iterative scheme is no exception. However, 
our method does provide a formal mean to quantitatively 
explore this aspect of the analysis, which was somewhat 
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overlooked in several previous works. We shall further 
refer to this issue later in this paper. 

The expression we obtained above for Vl/J [Eq. (52)] 



corresponds to the case of an initial pulse which is com- 
pact with respect to retarded time u. This pulse repre- 
sents a wide family of physically reasonable initial data. 
Yet another realistic initial setup involves the presence 
of a static field outside some radius (say, the surface of a 
spherically collapsing object) up to some moment of time 
(say, the onset of collapse) . The corresponding late time 
behavior of at null infinity can be easily inferred from 
our previous results, as follows. 

The static flat-spacetime solution to the Klein-Gordon 
equation (|l7|), regular at infinity, reads vpstat _ ^ r - J ; 
where fj, is a constant representing the strength of the ini- 
tial static field. The corresponding go function, namely 
the solution to Eq. @ with T(u) = ^ tat (v = 0), is 



2 l [l\/(2l)\]fi. In order to calculate \&i in this case, 



s-. St £L"t 

9o 

we need only to substitute for g in Eq. (|47|), while chang- 
ing the range of integration to f™^ du', u\ being some 
retarded time after the initial pulse ceases on the initial 
ingoing ray. Doing so, we find, to the leading order in 
M/u and in u\/u, 



Si°(«)[stat] 



>i+i 



M[l\/(2l)\]{iu 



which leads to 



¥J°(u)[stat] = (-2) 



r\ 2 



(21) 



■Mfiu 



-i-i 



(53) 



(54) 



This result is again identical with the one derived by GPP 
(Eqs. (12) and (18) in j|) for static initial data. 



VII. HIGHER TERMS OF THE EXPANSION 



(with t> N -i(v' < 0) = for any N > 1). 

With the explicit form of go and SV, and following 
the same integration-by-parts procedure as in deriving 
Eq. (p7|), the last equation can be put into the form 



9n(u,v) 



i 



Z < I I (ni _ q,l\2l — k+3 



x [a[ + b{ ln(v' - u')) g N -i(u', v>) B[ri - r.(JZ)], (58) 

in which a l k and b l k are constant coefficients. 

We observe that the analysis of for N > 2 in- 
volves two additional technical difficulties, which were 
not present in the calculation of First, the source 
wave is no longer compact with respect to retarded time, 
but rather it extends infinitely to the future. This means 
that ifftf will be influenced by sources located at any 
value of r, not only at null infinity. (However, as we 
discuss later, the small-r sources seem not to affect the 
asymptotic late time behavior at null infinity). The 
second technical difficulty is the fact that for N > 2, 
the source functions gpf-i(u', v') are also u'-dependent. 
Thus, in order to calculate f° r N > 2, one must first 
be provided with the form of gN-i at o,ny value of v, not 
merely at null infinity. 

To proceed, we shall first refer specifically to the case 
N = 2, in order to demonstrate that the dominant con- 
tribution to at late time comes from ^n—i sources at 
null infinity ; Sources at small r- values do not affect the 
behavior of the wave at null infinity to the leading order 
in M/u. Under this assumption, we shall than calculate 
this dominant contribution, to obtain the late time de- 
cay pattern of each of the functions 4" at at null infinity. 
Numerical support shall be presented. 



We now turn to evaluate the higher terms (N > 2) of 
the iterative expansion, each defined by Eqs. (O — fl6| ). 
These terms may be formally constructed in an inductive 
way by means of the 'flat space' Green's function derived 
in sec. [v|: 

^n(u, v) = 

du' / dv'G(u,v;u',v')6V(u',v')y N -i(u',v'). (55) 



Manifestly, this solution satisfies the initial and boundary 
conditions specified in Eqs. fll5| ) and (fl6|), respectively. 

In analogy to the treatment of [Eqs. ((39|) , (f4(i|)l , we 
can write 



V N (u,v) = — { [g N (u,v)\ 



(56) 



in which 



gN(u, v) = - du' dv'g G (u; u', v')8V(u', v')^ N -i(u', v') 

J uo J u 

(57) 



calculation of ^ 

To analyze Eq. ( |5q ) for N — 2, we separate the do- 
main of integration into four regions, labeled I — IV, as 
indicated in figure ^|. In what follows we show that the 
dominant contribution to ^^(u 3> M) comes only from 
sources in region I, which is a distant (large r value) re- 
gion of spacetime. We shall thus call this zone 'the main' 
region of integration. An upper bound will be set on the 
contribution from each of the other regions, II — IV, to 
verify their relative negligibility. 
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]U =u 



II a III 



u'=u/2 



u'=(Mu) h 



IV 

I 



1 u =u 

FIG. 6. Range of integration for VE 1 ^ (u). The contribution from 
each of the regions I — IV, indicated in this sketch, is evaluated in 
the text. 

First, consider region II, which is covered by u < v' < 
2u and u/2 < u' < u, with the sphere r' < R excluded. 
In this region of integration, the following upper bound 
is applicable to the source function gi(u',v'): 



(y> _ U 'V + 1 

\ 9l (u',v')\<CM [V , > ln(v'-uo), 



(59) 



where C is a positive constant. To prove this, set N = 1 
in Eq. (p8[), with the replacements (u,v) — > (u',v') and 
(vf, v') — > (it", v"), and with the u" integration cut off at 
u" = [since go(u" > 0) = 0]. Then use the fact that for 
any < k < I, 



(v" 



u 



ll\l—k 



(v" - u") 



//\2l-k+3 



< 



jv" - u') 
(«' - u") l+3 



< 



(v" 



u 



(«') 



'V+3 



(60) 



where the first inequality arises since v" > vl , and the 
second is due to the fact that u" < 0. The inequality 
( p9| ) then follows from (v" — u") < (v' — uq) and after 
integrating over v" . 

Now, to set an upper bound to the contribution from 
region II to ^f, consider Eq. (|58|) for N = 2, with the 
double integral replaced by f™, 2 du' J u u dv' (see figure ||) . 
Using the inequality (|59|), together with (u — u') < (v' — 
u') and (v' — u) < (v' — u'), one finds that 



|02?/(«)I < 

CiM 2 ln 2 u 
C 2 M 2 u- 3 In 2 



du' 



u/2 



2M dv >y- u :r < 



m ax { u , u ' + 2 R } 



'V+3 



(61) 



in which g^jj(u) stands for the contribution from region 
II to <?2°, an d where C\ and Ci are positive constants. 

Next, consider region III. The contribution g%ju(u) 
arises from the integration § u / 2 du' J„ dv' . The upper 



bound we have set above on gi(u',v') , Eq. (59), in no 



longer efficient at large r values. Here, we use a second 
upper bound, which may be is easily derived from Eq. 



\gi(u',v')\ <C 3 M(u')- 2 lnv', 



(62) 



where C3 is yet another positive constant. Provided with 
this upper bound, we observe that 



\9nn{u)\ <dM 2 / du' 

J u/2 



1 1 2 -/ 

A ' 111 V 

2u aV (u') 2 (v'-u'f 



< 



< C 5 M 2 U - 3 In 2 u, (63) 



with C4 and C5 being positive constants. 

We proceed by evaluating the contribution <?£fy(it) 
from region IV. Using again the inequality (|62|), we find 
that 



|ff§?v(«)| < C 6 M 2 



u/2 

du' 

'uM 



dl' 



In 2 v' 



< 



(u') 2 {v' -u' f 
< C 7 M 2 u~ 5/2 In 2 u, (64) 



where Cg and C7 are still another positive constants. 

Finally, we turn to calculate the contribution to (u) 
due to sources in region I, that is at u < u' < y/Mu. 
To that end we write the source function g\ in the form 



gi (u',v')=g<?(u')+Ai(u',v'), 



(65) 



in which Ai(u',v') vanishes at null-infinity. In appendix 
^ we show that the contribution to due to this in- 
dependent part of gi is bounded from above at large re- 
tarded time by 



|02aWI <C 8 M 2 u- 3 \n 2 u 



(66) 



where Cs is a positive constant. 

It remains to calculate the contribution to g^i due to 
g^°(u). With the source wave being a function of retarded 
time only, the analysis is then exactly the same as the 
above analysis for ^J . In analogy to Eq. (^), we thus 
immediately obtain, to the leading order in M/u and in 
u /u, 



g^ = -2Mh u 
in which we have defined 

fllu 



g?{u')du' 



g?{u')du' 



(67) 



(68) 



The first equality, accurate to the leading order in M/u, 
stems from the fact that the integrand is convergent at 



00 (since gi 



at large u) 



Collecting the above results [Eqs. (pip, (p3), (f 
and $57])] we conclude that 



nee 



• The dominant contribution to <?£°(u 3> M) [hence 
also to ^g°(u ^> M)} is only due to sources located 
in the 'main' region, namely at r > M; Sources at 
small r values have a negligible effect on <?£° at late 
time. 
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• To calculate the dominant contribution to g%° (u 3> 
M), one is allowed to replace gi(u', v') in Eq. (|5g[ ) 
by gf(u). This does not affect to the leading 
order in M/u. 

Therefore, Eq. ( p7| ) approximates the 'overall' late time 
behavior of 52 at null infinity. For the scalar wave itself 
we shall thus have, by Eq. (56), 



^f{u > to) = 2(-l)' +1 (Z + l)\Mh u 



-1-2 



(69) 



which is analogous to Eq. (p2Q. 

To summarize, we have found that \&2° poses the same 
late time behavior as namely a u~ l ~ 2 inverse power- 
law decay (for the compact initial data setup). To find 
out the relative amplitudes of these two terms at late 
time, one must calculate the ratio of the ir coefficients, 
ii/io- This will be done in section VIII, after we first 
analyze the behavior of the general term of the iter- 
ative expansion. 



The N'th term of the iterative expansion 

We would now like to understand the late time be- 
havior of each of the higher order terms (N > 3) of the 
iterative expansion. In principle, this may be done in an 
inductive way, using Eq. (|58|). We shall proceed by as- 
suming that, as was demonstrated above for ^ and T f 2 , 
the dominant contribution to at late time (for any 
N > 1) is only due to sources at null-infinity Phrased 
differently, we adopt the assumption that late time radi- 
ation at null-infinity originates predominantly from scat- 
tering off spacetime curvature at null-infinity. It means, 
in particular, that it is only the asymptotically-far re- 
gion of spacetime whose structure affects the asymptotic 
late-time radiation at null infinity. 

Accordingly, in order to obtain the late time form of 
^jv at null- infinity, we will use Eq. (58) with the source 
function gpf-i(u', v') replaced by g'^_ l (u'), and with the 
upper limit of the ^'-integration set to vf = yuM. This 
will presumably provide the correct form of ^^(u ^> M) 
to the leading order in M/u, as was explicitly indicated 
in the case N — 2. The treatment of the N > 3 cases is 
then fully analogous to that of ^2, and the generalization 
of Eq. (|6^) is straightforward: 



M) 

in which we define 
Ijv-i = 



-2M/jv-i u 



(70) 



(71) 



and which is accurate to the leading order in M/u and 
in uq/u. By Eq. (Jog) this finally leads to 



V%(u > m) = 2(-l) (+1 (Z + l)!MIjv-i u 



-1-2 



(72) 



We conclude that each of the terms ^n^i in the iter- 
ative series has a similar late time behavior at null infin- 
ity, that is a u~ l ~ 2 inverse-power law decay. Numerical 
experiments firmly support this observation, as demon- 
strated in figure ^. Consequently, one finds that, in gen- 
eral, to obtain the correct overall amplitude of the scalar 
wave at null infinity (even merely to the leading order in 
M/u), one must sum all contributions from the various 
terms ^n>i- If this sum turns out to converge — then the 
overall power-law should be u~ l ~ 2 . 



^1=1 



u/M 




-.1=0 

\ 
\1=1 



10 10 10 10 

u/M 




which generalizes Eqs. (p2j) and (|69| 



u/M 

FIG. 7. Late time tails of the iterative expansion terms. Pre- 
sented on a log-log scale are numerical results obtained for 
V&2 and >l/3 at v = 80000Af (approximating null infinity), for the 
I = 0,1,2 modes. Compact initial data for the propagation has 
been specified between u = —5M and u = — 10M, and the radius 
of the shell has been set to R = 3M. The results demonstrate the 
u~ l ~ 2 late time decay rate predicted by the analytic calculation. 

A considerable simplification is achieved when the case 
|uo| 3> M is considered. This corresponds to initial setup 
in which the pulse is specified where spacetime is ap- 
proximately flat (far outside the highly curved region). 
In this case, analytic arguments as well as numerical ex- 
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periments, both to be presented below, suggest that (i) 
the iterative series converges at null infinity at late time, 
and that (ii) the dominating term of the expansion there 
is if>x, which well approximates the "complete" wave \F 
(The contribution from the rest of the terms is smaller 
by order M/uq). 



VIII. CONVERGENCE OF THE ITERATIVE 
EXPANSION 

To discuss the convergence properties of the iterative 
expansion, we shall concentrate on the case |uq| 3> M. In 
what follows we argue that the ratio of any two successive 
terms of the iterative expansion, when evaluated at null 
infinity at«> M, is given to the leading order in M/uo 
by 



N-\ 



- oc Minor 1 In 



wo 



(73) 



JV 



(for N > 1) with a proportion factor depending on N. 
Below we present numerical results in a firm support of 
this relation. But first, we show that the validity of the 
ratio (|7^) is also suggested using analytic considerations. 

Throughout the following discussion we shall refer to 
the limit R/u a — > for simplicity. To that end we need to 
assume that the limit ^^(R/uq — * 0) exist, and that the 
functions *F^ are each independent of R to the leading 
order in M/u and in M/uq (this is supported by numer- 
ical examination). 

To evaluate the ratio 'FJ/ , ^/^^ we shall have to figure 
out how the various functions \P n scale with respect to 
the parameter |tto|. To that end we first show that each 
of the functions gjs;{u,v) can be written in the form 



g N (u,v)=u Q (ln|«o|) fN(u,v), 



(74) 



where /at are certain functions of only the dimcnsionlcss 
variables u = u/uq and v = v/uq. This form, which is 
valid globally (for any values of u and v) to the lead- 
ing order in M/uo, may be verified using mathematical 
induction: By Eq. ( p^ ) we learn that Eq. j7^ ) holds for 
go-*** Following the inductive procedure, we now assume 
that ([Fi]) is valid for a certain value of N. We then obtain 
for gN+i [using Eq. ( |58|) with x = u'/uq and y = v'/uo], 

g N+ i{u,v) = 



***To explore the scaling of go with respect to uq, it is con- 
venient to eliminate any scale parameters which may be char- 
acteristic of the initial data function T(u), by referring specif- 
ically to the case T(m) = 8(u — uo). Provided that any possi- 
ble scale characteristic of the initial date is much smaller then 
|uo|, this cannot limit the generality of our discussion. 



l-2-N 



k=0 



(\n\u \) N J2 / dx dy 



{y — u) l {u — x) 1 k 



t(0,«) 



(y - x) 



2i-fc+3 



x {a k + b k \n[u (y - x)}} f N (x,y). 



(75) 



For I Mo I 3> M this may be written (neglecting terms 
which are smaller by ~ In \ uq\) as 



g N+1 (u,v) = u l Q - 2 - N (ln\u \) N+1 fN + i(u,v), (76) 
where we have defined 



f N +i(u,v) = ^6 fe 



fe=0 



d.r 



dy 

max(0,u) 



(u - y) l (u - x) 
(y - ;r)2Z-fc+3 



l-k 



(77) 



Thus, by induction, Ec 

In virtue of Eqs 
N > 1) 



<Xli) holds. 



and (70) we then find that (for 



(78) 



which by Eq. (|7|) describes the scaling of the functions 
^^(u 3> M) with respect to uq. The ratio ( |7^ ) then 
follows immediately. [We comment that for the I = 
mode, no logarithmic factor is expected to occur neither 
in Eq. Q nor in Eq. @]. 

Eq. ( |73| ) suggests that the late time behavior of the 
scalar wave at null infinity is dominated by 'FJ , provided 
that I wo I 3> M. More accurately, it indicates that to the 
leading order in M/u, in Uq/u and in M/uq we have 



(79) 



The significance of this result stems from the fact that 
we have a simple analytic expression for 'FJ at late time. 

Numerical investigation, regarding the behavior of the 
functions 'F^ as related to the value of uq, provides a 
firm (though qualitative) support to both Eqs. ([73]) and 
([79|). Examples of these numerical results are presented 
in figures ||, ||, an d |l0| . 

Note that Eq. (|73|) by itself does not tell us whether 
the sum of the iterative expansion is convergent, as we 
do not know the form of the (iV-dependent) proportion 
factor appearing in this expression. However, the numer- 
ical experiments strongly suggest that the expansion does 
indeed converge, provided that |uo|/M is large enough. 
Also, the rate of convergence seems to increase as we take 
|uo|/M to be larger. These features are also apparent in 
figures H through [l0| 
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FIG. 8. Relative amplitudes of the functions \I/jv a t null infinity, 
for various values of the parameter uq . Shown are numerical results 
describing the 1 = mode of ^2 and ^3 at v = 10 5 M (approx- 
imating null infinity), for the cases m = -20M, -50 A/, —200 A/. 
Also shown for reference (in dashed line) is the 'complete' wave 
obtained by a direct numerical solution of the wave equation (Q). 
The results support the observation that the ratio |^l/]v/'I'jv+ij 
decreases as \uq\/M is chosen larger, and that for \uq\/M large 
enough, \E>J° becomes the dominant component of the 'complete' 
scalar wave if, as suggested by analytic considerations. Note the 
independence of the amplitude of 'I'i in the value of uq at late 
time, which agrees with the scaling rule ( fFsp for the monopole 
mode (where no logarithmic factor is involved). 




FIG. 9. The same as in figure g, this time for the I = 1 mode of 
the scalar field. The steep decreases in the amplitude of ^3 occur 
where this function changes its sign (see the comment in the text). 
In the case uq = — 200M, a "remnant" of such a sign-changing at 
u ~ lOOOAf "postpones" the development of the power-law tail to 
later times. 
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FIG. 10. Indications for convergence of the iterative scheme. 
Presented (on a linear scale) are the ratios (a) Wi/^ 1 , (b) — ^ — - , 
and (c) ^1+^2+^3 f or ; = o and for 1 = 1. Other parameters 
are R = 5M, uq = — 200M and v = 10° M (approximating null 
infinity). The results demonstrate the convergence of the iterative 
expansion for large \uq\/M values at large retarded time u. 



The amplitude of the various functions 'J? n may change 
its sign during the early stage of evolution (as apparent 
in figures |?] and ^J) . Our numerical experiments suggest 
that this kind of sign-changing occur at larger and larger 
values of u as TV and Uo increase, "delaying" the forma- 
tion of power law tails to later and later times. We could 
not rule out the (somewhat bothering) possibility that 
for any value of u and uq, there will exist some Nq such 
that for any N > Nq sign-changing would occur at re- 
tarded time greater than u. In that case, our analytic 
considerations [especially the ones leading to Eq. (73)] 
will not apply for all TV. Nevertheless, numerical exami- 
nation shows that the convergence itself is not disturbed 
by such sign-changing. 

Finally, a word of caution is in place: In our analytic 
calculations we were taking into account only the con- 
tributions to the functions * n due to scattering off the 



V~ In j-*' potential, which is the asymptotic form of SV 
at large distance. We have argued previously that con- 
tributions due to the other terms in the M/r* expansion 
of SV do not affect the behavior at null infinity at large 



u. This, indeed, is supported numerically, as demon- 
strated in figure [fl]. However, such terms (each giving 
an additive contribution to each of the functions 'Fjv) do 
affect the exact amplitude of the leading-order tail: For 
example, these terms give an additional contribution to 
the amplitude of WSf 3 , which (in terms of the M/uq ex- 
pansion) is of the same order as the amplitude of ^3°. 
Therefore, to examine whether the sum of the iterative 
series converges to the actual 'complete' wave (at 
u ^> M), one must take into account the 'complete' po- 
tential SV. (This, indeed, was done in obtaining the 
results presented in figures || through [To| ). 
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FIG. 11. Contribution from higher order terms in SV to VPS? * s 
negligible at u > M. This is demonstrated here for the 1 = 1 mode 
in the case R = 3M and uo = — 200M. Plotted are the ratios (a) 
VPi/vPj, (b) *2/*2' and ( c ) *3/*3> where the functions <S> N are 
those defined with the 'complete' potential SV, and the functions 
VPt, are defined with SV taken to the leading order in M/r*. 



IX. SUMMARY AND DISCUSSION 

The problem of scalar wave evolution in a curved 
spherically-symmetric spacetime can be formalized math- 
ematically as a two-dimensional initial-value problem, 
with a non-trivial effective potential. In this paper we 
incorporated the simple shell model in order to demon- 
strate the applicability, and explore the features, of a 
new calculation scheme to handle this problem. This 
new scheme is based on a special perturbative decom- 
position of the wave ("the iterative expansion") which, 
if fact, converts the original homogeneous non-trivial 
initial-value problem into an infinite hierarchy of inho- 
mogeneous initial-value problems, each having a simple 
Minkowski-like effective potential. The resulting initial- 
value problems can then be treated analytically in an in- 
ductive manner, using the Green's function in Minkowski 
spacetime. 

Unfortunately, we could not give a full analytic 
treatment concerning the convergence properties of our 
scheme. However, some analytic considerations, firmly 
supported by numerical analysis, indicated that the role 



1G 



of the "small parameter" of the iterative expansion is 
played by M/uq (although this parameter is not self- 
manifested in the wave equations). It was strongly sug- 
gested that the iterative expansion does converge at null 
infinity, provided that the support of the initial data for 
the wave evolution is confined to large distances (small 
M/|uo|). Moreover, our analysis indicated that under 
the above condition, the late time behavior of the 'com- 
plete' wave is well approximated at null infinity by 
namely the first element of the expansion manifested at 
late time. The last result is of special significance, since 
the late time form of \&i at null infinity can be derived 
easily, as shown in section [V]|. Hence, for |i*o| ^S> m our 
iterative expansion proves to be a simple effective tool 
for analyzing the late time behavior of the scalar field at 
null infinity. 

Here, it is important to distinguish between the issue 
of efficiency of the scheme to that of its usefulness. The 
iterative scheme can still prove useful even when M/|u | 
is too large for the expansion to converge efficiently (that 
is, for to be well approximated by merely ^i): In all 
cases when the scheme converges, it always predicts (by 
Eq. (fr^)) that the overall late time tail of has the 
form u~ l ~ 2 (for the compact initial set-up). The conver- 
gence may be 'efficient' (when |uo|/M is large enough), 
allowing one to obtain a simple analytic approximation 
for *S>°° by calculating Wi; or else it may be 'inefficient', 
in which case our scheme still provides a simple formal 
way to calculate fy 00 , by successive applications of the 
formula j55| ) (e. g. using numerical methods). 

We think the following is a reasonable physical ex- 
planation to the formal relation discovered between the 
smallness of the parameter M/|uo| (describing the dis- 
tance at which the initial perturbation originates), and 
the effectiveness of the iterative expansion. As already 
argued by Price 0], spacetime curvature at small dis- 
tances acts as an (almost) impenetrable potential barrier 
for the long waves of the scalar field. These long waves 
are the ones dominating the late time behavior of the 
field (the effective frequency goes to zero as t goes to in- 
finity). Since long waves coming from large distance can- 
not reach the highly-curved region (these waves are effec- 
tively scattered off the potential tail at large distances), 
their evolution takes place in the large-r region of space- 
time. Thus, for a perturbation originating at large dis- 
tance (|w | 3> M), late time behavior of the field should 
depend mainly on the structure of spacetime at large dis- 
tance; details of spacetime curvature at the highly-curved 
(small r) region should be less relevant for this behavior. 
This may no longer be true if the initial perturbation 
is specified at small r, thus allowed to evolve through 
the highly-curved region. Now, the iterative scheme has 
the mathematical structure of a perturbative expansion, 
with convergence properties that depend on how much 
can SV be considered a "small" amount. This function 
gets arbitrarily small as r gets arbitrarily large, but at 
small distances it is finite, and not small. However, the 
above reasoning implies that the form of SV at small dis- 
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tance seems to be irrelevant to the wave behavior at late 
time, provided that the initial perturbation originates at 
large distance (|uo| 3> M). The above argument supplies 
a physical explanation to our formal result that the ef- 
fectiveness of the iterative scheme increases as \uq\/M is 
taken larger. 

In this paper we have focused on analyzing the wave at 
null infinity. The results we obtained have a clear physi- 
cal significance with respect to a distant static observer: 
Along world-lines of constant r we have Au = At, where 
Am and At are, respectively, the retarded time and the 
(flat-space) static observe's time elapsed since the 'main 
pulse' of radiation had reached the observer. Now, for 
r very large, 'null-infinity' is in fact the region At <C r 
Thus, in accordance with our results at null infin- 
at times M < At « r such an observer observes a 
* cx (Aty l ~ 2 tail [or a * cx (At)~'~ 2 tail in the case of 
static initial setup]. 

There is yet another important outcome from the anal- 
ysis at null infinity: In the accompanying paper we apply 
a simple method to determine the late time behavior of 
the scalar wave outside the Schwarzschild black hole at 
any constant radius, down to the event horizon. This 
is made possible only if one is provided with the form 
of the wave at null infinity which shall be calculated in 
Schwarzschild using the same iterative scheme developed 
in the present paper. Thus, even thought the iterative 
scheme seems, in general, not effective at small distances, 
it still constitute, in our approach, a crucial step in es- 
tablishing a complete description of the wave behavior at 
late time. 

Generally speaking, we may point out the following 
advantages of the new technique: 

• Most important, the new scheme is the first one 
directly extensible to realistic rotating black holes, 
as shall be described in [ fl6[ 

• The iterative scheme provides a convenient formal 
framework for studying the effect of the non-trivial 
curvature-induced part ( U SV" ) of the effective po- 
tential on the evolution of the fields. This effect is 
isolated in our approach by explicitly splitting the 
wave into its 'trivial' (flat space) part V&o, and a 
curvature-induced part \& — ^o, in which each term 
^ n arises from N scattering events off the curva- 
ture potential SV. 

• We already discussed the fact that the iterative se- 
ries can be formally viewed as an expansion of the 
scalar field in the parameter M/uq. This feature 
provides a convenient formal way to explore and 
quantify the dependence of the wave behavior on 



t^The spacetime region At << r, for r very large, may be 
properly named (after Leaver ||) 'the astrophysical zone' of 
the waves. 
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the initial distance of the (compact) perturbation, 
(an issue somewhat overlooked in previous works.) 

Apparently, the iterative scheme developed in this pa- 
per can be directly applied to any asymptotically-flat 
spherically-symmetric spacetime (possessing an effective 
potential barrier for the waves). That would basically 
require only the modification of the curvature potential 
SV to that characterizing the new spacetime under con- 
sideration. The results of our analysis were consistent 
with the general assumption that the small distance de- 
tails of spacetime structure do not affect the behavior of 
the wave at asymptotically late time. (This was demon- 
strated for and if? 2 by an explicit calculation). Thus 
it seems reasonable to assume that the characteristics of 
this behavior would be completely determined by merely 
the large-r asymptotic form of SV. 

In the accompanying paper we use the iterative scheme 
to analyze the behavior of scalar waves at null infinity in 
the 'complete' Schwarzschild manifold. In this case the 
details of the analysis become somewhat complicated by 
the presence of the highly curved region near the horizon. 
Yet, we shall find that at late time the wave behaves 
essentially the same as in the shell model, in consistency 
with the above indications. 
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APPENDIX A: UNIQUENESS OF THE 
BOUNDARY CONDITION PROBLEM 

The purpose of this appendix is to show that the initial 
value problem for the scalar wave, with a boundary con- 
dition, as defined in Eqs. (^), (||) and (|l0|), has a unique 
solution. (The formal issue of uniqueness is not trivial 
because of the divergence of the potential function at the 
center of symmetry in the 1+1 representation.) It will 
also be verified that each of the components ^jv of our 
iterative expansion is uniquely defined. 



Uniqueness of ^0 

First, it will be shown that the solution for \I>o (the 
"Minkowski-like" first component of the iterative expan- 
sion), stated in Eq. (|l9|), is unique. To that end we exam- 



ine the most general solution to the wave equation (17) 
given by 



(» 
1 9o 



(v 



(Al) 



in which the coefficients A l n are those given in Eq. ( |2C 
and where g(u) and h(v) are (yet) arbitrary functions. 
(In flat spacetime, these two functions describe the out- 
going and ingoing components of the wave, respectively). 
It will now be shown that although the functions g(u) 
and h(v) are not uniquely determined by the initial and 
boundary conditions, the wave ^0 *s unique. 

Consider first the evolution of the wave between re- 
tarded times u = uq and u = (see figure ^| for refer- 
ence). In this region the wave equation possesses no ir- 
regular points, and thus existence and uniqueness of the 
solution ^0 is guaranteed relying on standard theorems 
(see, for example, p9|). Therefore, the solution stated in 
Eq. (19a) is unique at u < 0. 

We thus focus on the region u > 0, where we have 
to show that with the conditions ^q{u — 0) = and 
"f?o(r = 0) imposed on the general solution (Al), the 
trivial solution $?q(u > 0) = is unique. (Due to the 
linearity of the wave equation, to verify uniqueness it is 
enough to consider the case in which the wave vanishes 
on the initial ray u — OPL Thus the use of the condition 
^o( M = 0) = in this context is regardless of the fact 
that the wave actually vanishes along the outgoing 
ray u = 0.) 

Imposi ng t he condition $?q(u = 0) = on the general 
solution (Al) leads to a linear inhomogeneous ordinary 
equation for the function h(v ), analogous to Eq. (pl|). By 
analogy to Eq. ( p2|) , the general solution to this equation 
reads 



1 



1 ^^^j+iTt^T 1 
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' 7). 



(A2) 



n=l+l 



where c„ are I arbitrary parameters. Upon calculating 
the integral, we find that h(v) is just a polynomial of 
order 2L 



11 ^ 



(A3) 



with its first I + 1 coefficients given by 

c« = -5o n) (0) for0<n<L (A4) 

Next, we require that would vanish at 7% = (that 
is for u — v). By Eq. (Al) this requirement makes the 



'"For, if there exist two solutions at u > 0, both subject 
to the same initial and boundary conditions, then their dif- 
ference admits the same wave equation, with zero initial and 
boundary conditions. Then, showing that the difference is 
identically zero implies that the two original solutions must 
coincide. 
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functional equality g(u) — —h(u) com pulsory. Note that 
this result accommodates with Eqs. (A2,A4), allowing 
both conditions iff (u = 0) = and * (7* = 0) = to 
be satisfied simultaneously without forcing the functions 
g(u) and h(v) to vanish identically. 
We conclude that \&o has the form 



= E4 



Vl n) (u) - (-l) n V^>(v) 



(n). 



n=0 



(v — u) 



l — n 



(A5) 



where Vi are polynomials of order I [the (1+ l)th through 
(2Z)th powers of the polynomials vanish independently in 
the summation, in view of Eq. (p3|)]. Notice that formally, 
iffo vanishes along the ray u = for any polynomial Vi 
whatsoever, as the polynomial coefficients c„ are arbi- 
trary in our construction. 

The last remark is important, since it directly implies 
that "Jo vanishes anywhere at u > 0. For, given a specific 
value of retarded time u > 0, we may transform to the 
new coordinates x = u — u and y = v — u, in terms of 
which we have 



l—n 



n=0 
I 

n=0 



(y — x) l ~ n 



(A6) 



where Vi is still another polynomial of order I. The last 
expression vanishes at x — (for any i>) regardless of 
the form of the polynomial, and thus we find that \&o 
vanishes also along the ray u = u. Since the value of u 
is arbitrary, we are led to the conclusion that iffo must 
vanish identically at u > 0. Therefore, the solution jl9| ) 
for is unique. 



Uniqueness of the functions ty? 



It can be shown that no solution exist for the Green's 
function, which is continuous along u — v' (with the sin- 
gle point u = v = v' excluded) and, at the same time, 
satisfies the boundary condition. For, if we assume that 
G(u — ► (v') + ) = 1, then we shall have the function he 
admitting the form 



21 

h G {v) = Y j -,c n {v-v') n 



(A7) 



n=0 



[analogous to Eq. (A3)], where this time the coefficients 
Cn shall read 



-9g\v') + 1 



,{)<n<l-\ 
,n = I 



(A8) 



This cannot be accommodated with the requirement that 
^g(u) = —go(u), necessary for the boundary condition 
to hold. Thus the Green's function cannot be continuous 
along u = v' . 

Now, the Green's function can still obey the field equa- 
tion at u = v', provided that the discontinuity along this 
ray is constant in magnitude for all v > v'. This can be 
verified by writing G = G X 6(u — v') + G X 6{v' — u), and 
substituting into the equation for the Green's function. 
One thereby finds that for this equation to be satisfied 
at u = v', we must have G, v (u — > (v') + ) = 0, and thus 
G(u — > (v') + ) — k, where k is constant. By analogy to 



Eq. flA8D we find that in this case h^(v') = -g£ } (v')+k, 
which implies that the boundary condition is violated, 
unless k — 0. Hence, the Green's function must "jump" 
to zero right after retarded time u — v'. 

Then, by the discussion regarding w ^ conclude that 
the only solution to the Green's function in region II, sub- 
ject to the boundary condition at the origin and contin- 
uous everywhere except on the ray u = v', is the trivial 
solution G(u > v') = 0. 



,(0/ 



To show that the functions ^jvj for N > 1, are unique, 
it is enough to verify that the Green's function G is 
unique. Then, all function are unique by construc- 
tion. 

To discuss the uniqueness of the Green's function, we 
refer to a source point at (u',v'), as sketched in figure 
|3|. In region I (see the figure) existence and uniqueness 
are guaranteed by standard theorems ]l9[| . In r egio n II 
the Green's function admits the general form (Al), as 
do (both functions obey the same homogeneous field 
equation in this region) . In the expression corresponding 
to the Green's function we shall denote the two functions 
of u and v by gciu) and hc(v) respectively. 

The value of the Green's function along the outgoing 
ray u = v' , dictated by the evolution in region I, is unity. 
Hence, if the boundary condition G(r* = 0) = is to 
be satisfied, then there has to be a discontinuity in the 
Green's function at retarded time u = v = v' (the left 
vertex of the dark-colored rectangle in figure ||) . 



Uniqueness of ^ 

Now, consider the initial value (and boundary condi- 
tion) problem for the "complete" scalar wave iff in the 
shell model. To verify that a solution is unique, again 
we analyze the situation of zero initial data (on u = uq 
and v = 0), showing that in this case the null solution 
(iff = 0) is unique. 

At u < 0, existence and uniqueness of iff is assured, 
relying on standard theorems. Thus the trivial solution 
^(u < 0) = is unique. Particularly, we find that iff 
vanishes along the ray u — 0. 

Now, consider the region u > 0, v < 2R, represented 
in the diagram of figure |l2| by the triangle ABC. In this 
portion of spacetime (confined to the flat interior of the 
shell) iff obeys the same equation as if>o, and is subject to 
the boundary condition at = and to zero initial con- 
dition along the null segment AB (see the figure) . By the 
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above discussion (regarding ^o) we deduce that \& must 
vanish in the portion ABC. In particular, we find that 
VP is zero along the null segment BC. This, in turn, to- 
gether with the fact that the wave vanishes along the ray 
BD (sec figure), implies (again by standard theorems) 
that \& vanishes anywhere inside the region represented 
in figure |lj by the rectangle BDFC. Thereby we have 
shown that the wave vanishes anywhere between retarded 
times u = and u = 2R. 




FIG. 12. Uniqueness of the solution for *I/ (see explanation in 
the text). 

We can now proceed, referring next to the triangle 
CEG (see figure), to show that vanishes also between 
retarded times u = 2R and u = AR. The same procedure 
may be applied an arbitrary number of times, showing 
that the wave \I/ must vanish anywhere. Since the wave 
equation is linear, we have thereby shown that the solu- 
tion for 'J (with arbitrary initial data and zero boundary 
condition at the center of symmetry) is unique. 



n successive times with respect to v! . As explained pre- 
viously, compactness of go(u) assures that the resulting 
surface terms would all vanish. One is thus left with 



-Y 

n=0 



B l n (-l) n / du'g (u') 



gn 

du' n 



ln(v' - u') 
(u - u') 2 ~ n 



Summing first the n = and n = 1 terms yields 
b f° 9o{u') , 

tCLU , 



2(2l + l)J U0 {u-u'Y 



(B2) 



(B3) 



in which the vanishing of the logarithmic dependence 
arises from the equality of B l a and B[, which was also 
responsible for the vanishing of the contribution from 
scattering off the "r~ 3 " potential. 

The remaining terms in the sum in Eq. ( |B2| ) take the 
form 



I n— 2 



n=2 j=0 



l n\(n- 2)!(-l) n +J 



go("') 

j\(n -j- 2)!(n - j) J U0 (u - u')'- 



■du 



■du' (B4) 



(where the expression given for the sums over j and n 
is not too difficult to verify, using simple combinatorial 
manipulations) . 

The sum of two expressions given in ( |B3|) and (B4) 
corresponds to the logarithmic part of Eq. (|Bl|). We still 
have to consider the contribution of the non-logarithmic 
terms. Integrating these terms by parts n times with 
respect to u' (with the resulting surface terms again van- 
ishing), we find that only the n = and n = 1 terms 
survive, with their sum given by 



APPENDIX B: DERIVATION OF FOR I > 1 

In this appendix we give a detailed derivation of gf (u) 
[Eq. (|47j)] for the case I > 1. The starting point to our 
It has been shown that the r~ 3 



calculation is Eq. 
potential [namely the term proportional to a in Eq. (41)] 
gives no contribution to <?J° . We therefore set a = 0, 
then integrate Eq. (|4|) by parts once with respect to v'. 
Using Eq. ©, this yields 



ln(f/ -u') + J2 



'g^(u')(u-uT- 2 

n 



m—Q 



21- n-m + 2 



(Bl) 



with the coefficients B l n given in Eq. 

To proceed, consider first the logarithmic part of Eq. 



( Bl ) . We integrate by parts each of the terms in this part 
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; m=0 



1 



1 



21 



m + 3 21 — m 
'"' 9o(u') 



x 



ad 



(u — u') 2 



-2 
du', 



which is 



i(i + i)(2i + i)j U0 (u-u'y 



■du'. 



(B5) 



(B6) 



We finally obtain g^°(u) by collecting the expressions 
@, (H) and (H|), and setting b = 8Ml(l + 1) [by Eq. 
(m\. This produces Eq. @ 



APPENDIX C: 
CONTRIBUTION TO gf DUE TO Ai(w',w') 

The purpose of this appendix is to prove Eq. (|66|). This 
equation sets an upper bound on the contribution to gf 
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due to the ^-dependent part of g±, namely Ai(u, v) (de- 
fined in Eq. (pq)). 

We first show that the following upper bound is appli- 
cable to Ai (it, v) at large r: 



I Ai (it, v) | < CMrZ 2 In f, 



(CI) 



in which C is a positive constant. To verify the validity of 
this inequality, examine the explicit form of Ai, derived 
directly from Eq. (|58|), 



Ai (it, v) = / c?it' 

fc=0"'"° 



,{v' -u)\u-u') 



i\l — k 



(v' — u') 21 ^^ 3 
x [a{ + b{Hv'-u')}g a (u'). (C2) 



Now, since (1/ — it) < (1/ — it') and (it — it') < (1/ — it'), 
we find that there exist positive numbers C , C" , C" and 
C, such that 



|Ai(u,«)| < 

Cm / du' 



>u 

C"M I du 



efo — — \g (u )| < 



" , ln(i)-fi') 
(u - w') 2 



(1/ - it') 3 

ffSV)l < 



(l) — 1i) 2 



CMr" 2 Inn 



(C3) 



which confirms Eq. (CI). 

We are now in position to evaluate the contribution of 
Ax to g%°(u > M). By Eqs. @ and @ we deduce 
that this contribution (in absolute value) is bounded from 
above by 



fe=0 



(C4) 



Since (v' — u) < (u' — it') and (11 — u') < (1)' — it'), we find, 
in turn, that the last expression is bounded from above 
by 
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(1/ — u'y 



{u - u'Y 



< CM 2 u~ A In 2 ii, (C5) 



for 11 ^> M. [Here, C , C " an d C are some positive 
numbers, other then in Eq. (C3)]. 

The inequality (pq) is thereby proved. 



21 



